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Q> ' Abstract 

0> \ 

■ As is known from QED, a possible solution to the ghost-pole trouble can be 

$-h ! obtained by imposing the Q 2 -analyticity imperative. Here, the pole is compensated 

by the a non-analytic contribution that results in finite coupling renormalization. 
^ ■ We apply this idea to QCD and arrive at the Q 2 analytic a s (Q). This solution 

corresponds to perturbation expansion, obeys AF and, due to non-perturbative con- 
tribution, has a regular IR behaviour. It does not contain any adjustable parameter 
(N| ■ and has a finite IR limit a s (0) which depends only on group symmetry factors. In the 

one-loop approximation it is equal to 47r//?o — 1.40 and turns out to be surprisingly 
\ stable with respect to higher order corrections. On the other hand, the IR behaviour 

CO 1 of our new analytic solution agrees with recent global low energy experimental esti- 

mates of a s (Q 2 ). 

G\ ; 

1. We recall first some results obtained in QED about 40 years ago. The QED effective 
coupling a(Q 2 ) being proportional to the transverse dressed photon propagator amplitude 
is an analytic function in the cut complex Q 2 plane and satisfies the Kallen-Lehmann 
spectral representation. The "analytization procedure" elaborated in papers |], consists 
of three steps: 

•*-| ■ (1) Find an explicit expression for qirg{Q 2 ) in the space-like region by a standard RG 

improvement of perturbative input. Continue this expression to the time-like Q 2 domain. 

(2) Calculate the imaginary part of olrg{—Q 2 ) on the cut and define the spectral density 
Prg{°i ol) = Im a RG (-a, a). 

(3) Using the spectral representationwith p R Q in the integrand define an analytic 

OLan(Q)- 

For one-loop massless QED, this procedure produces || an explicit expression (see 
Eq. (2.6) in Ref. or analogous QCD Eq. (0) below) which has the following properties: 

(a) it has no ghost pole, 

(b) considered as a function of a in the vicinity of the point a = it has an essential 
singularity of the exp(— 3n/a) type, 

(c) in the vicinity of this singularity for real and positive a it admits an asymptotic 
expansion that coincides with usual perturbation theory, 

(d) it has a finite ultraviolet asymptotic limit, a(oo,a) = 37T, which does not depend 
on the experimentally input value a ~ 1/137. 

The same procedure in the two-loop massless QED approximation yielded a more 
complicated expression with the same essential features. 
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2. To use the same technique in QCD one has to make two observations. First, since 
a s (Q 2 ) has to be defined via a product of propagators and a vertex function, validity of the 
spectral representation is not obvious. However, this validity has been established in the 
Ref.|| on the basis of analytic properties of the vertex diagrams. Second, for QCD with an 
arbitrary covariant gauge running of the coupling and gauge parameter are interconnected. 
For simplicity, we assume that the MS scheme is used. 

In the following we use the spectral representation in the non-subtracted form 



The usual massless one-loop RG approximation yields the spectral function 

.(1) 



Substituting p { R Q into spectral representation Eq. ([[]) we get 



Uq 2 ) = l 



1 A 2 
+ 



lnQ 2 /A 2 A 2 - Q 2 



(2) 



where we have used the QCD scale parameter A. The "analytic" running coupling Eq. (|2D 
has no ghost pole and its limiting value 

a?>(0) = 47T/A, (3) 

does not depend on A being an universal constant which depends only on group factors. 

We have become accustomed to the idea that theory supplies us with a family of possible 
curves for a s (Q 2 ) and one has to choose the "physical one" of them by comparing with 
experiment. Here, in Eq. (H) the whole bunch of possible curves for d s (Q 2 ) corresponding 
to different A have the same limit at Q 2 = as shown on Fig. 1. For rif — 3 it is equal to 

5^(0) = 4tt/9~ 1.396. (4) 

Another feature of Eq. (Q) is the fact that its correct analytic properties in the IR 
domain are provided by a nonperturbative contribution^ like exp(— l/a/? ). 

To investigate the stability of the result, Eq. (f|), with respect to the next loop correc- 
tions, we have considered the two-loop approximation to a s (Q 2 ) in the form 

-( 2 ) 1 t i & h A a mo 38 



l RG 



and 



A) [L + hlnL] ' A 2 ' /?o ' 3 



R(£) = £ + 6i ln (V^ 2 + 7T 2 ) , = 7T + 6i arccos 



7T Z 



x The connection between "Q 2 analyticity" and "a nonperturbativity" has been discussed in the Ref. 
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The limiting two-loop value a^(0) is also specified only by group factors via /? and 
/3i . Surprisingly, this value found by numerical calculation practically coincides with the 
one-loop result. For the MS scheme in a three-loop approximation a@'(0) changes by a few 
percent. Thus, the value a s (0) is remarkably stable with respect to higher loop corrections 
and is practically independent of renormalization scheme. 

3. To fix A we use the reference point M T = 1.78 GeV with a s (M 2 ) = 0.33 ± 0.03 §. 
Corresponding solutions a^ =1,2 ^ (Q 2 ) are very close to each other for the interval of in- 
terest Q 2 < 10 GeV 2 . For instance, at Q 2 = 10 GeV 2 we have ^(lO) ~ c^(10) = 
0.267, (10) = 0.265 . Here, again we used rif = 3 as the average number of active 
quarks in the spectral density. This seems to be reasonable in the IR region. 

For a more realistic description of the evolution of a s (Q 2 ) in the Euclidean region 
3 GeV < Q < 100 GeV, one should take into account heavy quark thresholds. Using the 
explicitly mass-dependent RG formalism developed J7| in the 50's a "smooth matching" 
algorithm has been devised recently 0. This can be given to correct analytic properties 
while incorporating heavy quark thresholds. 

The idea that a s (Q 2 ) can be frozen at small momentum has been recently discussed 
in some papers (see, e.g., Ref. ||). There seems to be experimental evidence indicating 
behaviour of this type for the QCD coupling. As the appropriate object for comparison 
with our theoretical construction we use the average 



A(k) 



dQa s {Q 2 ,K). 



(6) 



"Experimental" estimates for this integral are A{2 GeV) = 0.52 ± 0.10 GeV |J and 
A{2 GeV) = 0.57±0.10 GeV HlQ| . Our one-loop results in case Eq. (H) for A are summarized 
below for a few reference values of a s (M 2 ). 



a s (M 2 T ) 


0.34 


0.36 


0.38 


A{2 GeV) 


0.50 


0.52 


0.55 



Note here that a nonperturbative contribution, like the second term in l.h.s. of Eq. (^|), 
reveals itself even at moderate Q values by "slowing down" the velocity of the a s (Q 2 ) 
evolution. For instance, in the vicinity of c and b quark thresholds at Q = 3 GeV it 
contributes about 4% which could be essential for the resolution of the "discrepancy" 
between DIS and Z data for a s (Q 2 ) . 

4. We have argued that a regular analytic behaviour for a s (Q 2 ) in the IR region could 
be provided by nonperturbative contributions which can be considered as a sum of powers 
ofA 2 /Q 2 . 

Probably, our most curious result is the stability of a "long-range intensity of strong 
interaction", a 8 (0) , as well as the a s (Q 2 ) IR behaviour that turns out to agree reasonably 
well with experimental estimates. 

On the other hand, the nonzero a s (Q) value evidently contradicts the confinement 
property. To satisfy this, one should have 

fia(0) =0 



as it has recently been emphasized by Nishijima [11 in the context of the connection 



between asymptotic freedom (AF) and color confinement (CC). 
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It is possible to correlate this type of the IR limiting behaviour with the RG-improved 
perturbative input and Q 2 -analyticity by inserting a Castillejo-Dalitz-Dyson zero into our 
solution (see Ref. |12|]). Such a generalysed analytic solution will contain additional pa- 



rameters. In this construction there is no evident relation between AF and CC. Here, CC 
is provided by nonperturbative contributions. 
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